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ABSTRACT 


A theoretical investigation is developed herein in order to 
estimate the permanent transverse deflections of rectangular plates 
when subjected to large dynamic loads. The influence of finite de- 
flections or geometry changes and of strain-rate sensitivity is retained 
in the analysis but elastic effects are disregarded. The particular 
case Of a fully clamped plate acted on by a uniformly distributed 
pressure pulse is studied in detail. It is observed that good agreement 
between the theoretical predictions and experimental results has been 
obtained for rectangular plates made from strain-rate sensitive 


material. 
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NOTATION 


B Semi-width of plate 

D Material constant defined by equation ( 2 ) 

H Thickness of plate 

I pot 

1/5 

I Oe al Hp Cc) 
1, Semi-length of plate 

SoH" 

Mo Static limiting plastic moment = a 
M'o Dynamic limiting plastic moment, defined by 


equations (4a-d) 

No Static limiting plastic axial force = 60H 

IN=O . Dynamic limiting plastic axial force, defined by 
equations (4a-d) 


Durations of parts of the response 


Vo Initial velocity 

Wm Permanent transverse deflection 

Wf Permanent transverse deflection if 2 Hal 
p Material constant defined by equation (2) 
pe Magnitude of static collapse pressure 


pi, p2 External pressures per unit surface area of a plate 

po Magnitude of dynamic pressure pulse 

t, to, tg Time 

Ww Transverse deflection of the midplane of a plate 

x,x',y Co-ordinates defined in figure 2 

Z Co-ordinate which is perpendicular to the mid-plane of 


an initially flat plate 


b B 


L 

E 1] Direct strain 
ee, 

ne iG 





sats 5 a tenes 





‘ ; ye o>s- 


(") 


fo 


Relative angular rotation rate across a line hinge 
Curvature 
Mass per unit area of plate 
tang 
Non-dimensional times 
Yield stress 
Duration of a rectangular pressure pulse 


angle defined in figure 2 
ot 


strain at middle axis of plate's section 
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INTRODUCTION 


In reference {1] Jones investigated the influence of finite 
deflections on the behavior of arbitrarily shaped plates and beams 
which are subjected to large dynamic loads. In his analysis Jones 
disregarded any elastic effects, strain-hardening and rate sensitivity. 
It has been shown that elastic effects can be neglected in the computa- 
tion of final deformation if the kinetic energy input is appreciably 
larger than the maximum elastic strain energy the plate can absorb. 
With respect to material properties, experiments have indicated that 
strain rate sensitivity is by far the most important property neglected 
in rigid- perfectly plastic methods. Thus, while experiments with rate 
insensitive materials, suchas Aluminium A16061T6, show remarkable 
agreement with theoreticai results, in the case of mild steel, a 
notoriously rate sensitive material, significant disagreement between 
theory and experiment takes place, the net result being an underestima- 
em on the ability of a plate to sustain a given dynamic load. Cowper 


3,8] and Symonds and Jones [4] have developed an 


and Symonds 
approximate theoretical method to examine the influence of rate sensi- © 
tivity in beams. Their results show drastic improvement, as compared 
with experimental evidence, over the non-rate sensitive methods. pe) 
A general approximate theoretical procedure, which retains both 
finite deflections and strain rate sensitivity, is presented herein for 
the dynamic behavior of arbitrarily shaped, rigid-perfectly plastic 
plates. This procedure is essentially an extension of the one developed 
by Jones [1] in order to include strain rate effects. The particular 


case of a fully clamped rectangular plate subjected to a uniformly 


distributed dynamic pressure pulse is studied in more detail. 





Pe PROXIMATE ANALYSIS 


In reference [1] Jones shows that if a plate (or beam) is 
divided into a number of rigid regions separated by ¥ straight line 


hinges, each of lengthlm, then an energy balance yields the following 


relation: y 
) (pe = ye) lec Ve Pa ee. M) Gancten (1) 
A m=1 lm 


where w is the transverse deflection at the hinge. 

The stresses at the hinges due to the axial force N and bending 
moment M act ona plane which is parallel to the hinge and transverse 
to the midplane of the plate. Om is the relative angular rotation rate 
across the hinge and is supposed to be infinitesimal. The quantity 
(Nw - M) 8m, on the right hand side of equation (1), represents the 
internal energy dissipation per unit length of a hinge. This dissipation 
function will depend on the boundary conditions of a plate and on the 
yield condition used. 

For a strain rate sensitive material Symonds and Jones in 
reference [3] represent the dependence of the dynamic lower yield 


stress 6 (E) on the plastic strain rate E by the following: 


Cu 
Se) 2 || Saas (2) 
So 7 D 


where 6% is the static yield stress, and p and D are material 
constants. This relation is totally empirical and approximate, pro- 
viding good representation of test results up to strain rates in the 
neighborhood of 1,000 sec7!. If formula (2) is used to derive dynamic 


plastic bending moment and axial force for rate sensitive behavior, 


[4] 


assuming the Navier hypothesis valid, Symonds and Jones show 


that: 
(see next page) 
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Mo'! a ie H2 (4a) 
N _ Z 
No! | : H (4b) 
_ 2 Ee alas 
where Mo [ 1 Opti ( oD )'/p ] Mo 
No' =[ 1+( ) “hp ] No 
and Z = 
Vi ee a si 
Ms — 2p+4 2Z 
> | 
o 
where aie) = fal + (22) Py Mo 
D 
P 
No = [1+(£)"] No 
ox H 
nod a > a 
K Seeulvatire rate 
Zz = distance between middle and neutral axes of the plate's 


cross-section. 

The moment Mo' and axial force No' are to be interpreted as 
the dynamic equivalents of the static values Mo, No. 

When equations (4a) and (4b) are combined they give the 


yield condition: 





\mr\ Ny 2 YIELD 
Mo! io CONDITION (5) 


- }]- 





This condition is the maximum normal stress yield criterion, para- 


bolic in shape, as is seen in figure 1. 





figure 1 
It should be pointed out that equation (1) has been derived for 


arbitrarily shaped plates and since it is an energy balance equation, it 
can be used with any kind of edge conditions, provided that plastic 
hinges are allowed to form around the boundary cf the plate, if 
necessary, as well as the interior. 

Ifa kinematically admissible collapse mechanism with straight 
line hinges is postulated, then equation (1), combined with the 
appropriate relations for (Nw - M) Om, can be solved to give an 
estimate of the influence of finite-deflections and strain rate on the 
deflection-time history of rigid, perfectly plastic plates (or beams) 


loaded dynamically. 


Peel Vic BanAVIOR OF A FULLY CLAMPED RECTANGULAR PLATE 


From infinitesimal plasticity analyses it is evident that the 
shape of the displacement field of a beam or plate for small dynamic 
loads is the same as the static collapse velocity profile [1,2]. In 


[2.5] it is seen that the permanent 


addition, from experimental results 

deformed profiles of rectangular plates loaded dynamically are similar 

to the shape of the velocity field used by Wood [6] in static analysis. 
Consider a rigid, perfectly plastic rectangular plate of length 


2L and width 2B, fully clamped around its outer boundary, as shown in 
figure 2. The plate is subject to a uniformly distributed dynamic load 


with the pressure-time history shown in figure 3. 
oe 





This external loading may be expressed in the form: 


} 

Do 4 ie fOr oXt KG ; (6a) 

Pa = 0 OO (G Kt << Aa (6b) 
where C and (T +T ) are the durations of the pressure pulse and 


Pieepiate response, respectively. 


pw 


fe 


wag 
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figure 2 fioure 3 
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Assuming Wood's displacement field, as discussed above, we 


laeiiaere 


Mee etanebie 5c") rn 
= [>= a (7a) 
B tang 
for region I 
and: 
eee ane cea y) (7) 
B 


for region II. Regions I and II are indicated in figure 2 andi = 1 
refers to the deflections during the time imterval 0 x t <G » while 
1 = 2 refers to Feat < eC. 

In order to use equation (1), its right hand side is examined 


mimst. thus: 
(see next page) 


wise 





i B 





COS 
Dein dim = 4 { Nearer an 
5 sing 
m=1 lm O 
L-B tane L B 
+2 - M)9m dx+4 | (Nw - M)@m dx +4 | (Nw - M) 9m dy ey 
O O 


The evaluation of the right hand side of equation (8) requires the use of 
equations (7a, and 7b) together with the yield condition and the 


relations: 


INV GS) Tee A No = 60H 


Values for the curvature rate K, rotation rate Om and hinge width to 


be used in each integral are given in table (1). 


Integral ¢ Integral Integral * Integral 
poet 
. sin 


maw oo! 
2Bl sin Bl tan 
ee 
Widt 


Table (1) 







It should be pointed out that in the third and fourth integrals w = Oe 
If the above information together with equations (6a), (7a) and 
(7b) are substituted into equation (1) then it is straightforward to show 


hat: 
poet (see next page) 


eA 





W 
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where: Q <K t 


2 - y 2 
+ h ta + = 
( 1 Ws hy) va p h,W, qd, (9) 











= Ms ee es hy = 2 , d, = 23? 
i or Pi By 
2 
__HBo (2-btand) | | _ B 
12 Mo (1+ b cote ) ae 
— [ 1 4 2 (1 - btan¢ )} 
Sil (1 +bcoth ) 
_ se . Pe -22Mo lt ecole 1) 
Po 7? oe (3-btane ) 
A 
ee 
Giant cote ) 
i A, b 
2 
omelets be cote: ) 
Crs C. si : ON oe @ nes 
_ 9 : 9 Sin 1° &3 ©0 
B 2 sin COS 
2 
Ss = _ 3C, tang + 3C)-C, afiees -C, - Cosa 
: B 2sin< cose 


Salas 





c. -2R— ,H_¥, 


oo 2 Bip 

_ 22 ae: 
ee 461 sine? © 
» 2p H yl 


3 2p+1 '2BlDtand’ P 


and taneb = \/3 = oe - b to ensure smallest upper bound to the 
actual static collapse pressure. i 
Equation (9) obviously is very difficult to solve analytically. 
Given the fact that the whole procedure herein is an approximate one, 
and aiming at providing the designers with relatively simple results, a 
drastic simplification will be made by setting p = 1, D = 500 sec - albeit 
Slesproper values of p and D for mild steel are 5 and 40 respectively. 
Thus, the highly non-linear equation (2) is linearized. Even so, 
however, it is hoped that the whole procedure will give sufficient insight 
into the behavior of strain rate sensitive plates. 
Now equation (9) may be recast as: 


Zz 
i 


+h) W+h,W,” _ (9a) 


+ 
Mi on on el 1 


A series solution to equation (9a) employing the method of undetermined 
coefficients, and which satisfies the initial conditions RS W, = O at 


t = 0, yields the results: 


and W =- W 


at the end of the first stage of motion (t = ), where: 








7 2 3. 3 44 
Ww, = d,@ Cie eee 6 2 Be Lb 
2 3 12 60 360 
d ht h.d.t?  h : fn does 
> Se Ara 2 © 4 Agh, 4,Tt 
60 420 2520 126 
re 5 
84 
2 2 3. 3 4 4 
ee et tet to 
1 1 2 6 24 120 
4 5 as 
— Aga, T e hed, © : ey 75 " 
20 120 720 
5 ES 
«, ee = eerie Se (9c) 
36 24 


A study of the second stage of motion ( TL a t £°T+T) 
proceeds in a manner similar to that outlined above for the first stage, 
but-with equations (9b) and (9c) as initial conditions, and with Y= 0. 


Thus the differential equation to be solved now is: 








e.6¢ ? a ) _ 
Wy + (h, W, gat) Wi + how, = dy (10) 
eet ee eed 

where d, a : ( ) = at 


Sige 





A series solution to equation (10) employing the method of undetermined 
coefficients and satisfying equations (9b) and (9c) gives the maximum 


value of the transverse deflection of the permanently deformed plate as: 


ee eee ce Ee EE 
a See SE ee 





H H H H H H H 
where: @, = °2 - 7 (hv +h wW.*)w, -1 nw.’ 
2 — 0 1" 17 — "34 
D 2 9 
cae, [&. (h. +h Fey eae a oe | 
3 7 De eh as) cs Pa aa | ieee 


4 =- [2 W.. W “ - +h W = 
4 [2a, 1 (3h, 1 h,) 3H, (hy h, Ww, at 


me ae y) 
Se 
oe, Wi) WY 


AG 
and i a i is the smallest non-negative root of the poly- 


nomial equation: 


3 3 2 2 4 
. + - 
Nome oe Pec pt Ww, 0 


4 (12) 


Equation (12) is obtained from the requirement that Wh Oat y= Ag 
where ty refers to second stage alone. 


Equation (11) remains valid provided Wm ¢ 1, because when 


S 
Wm oan 
Pie y 1 the dissipation expressions (Nw - M) 9m_ take on new 


values ( since then N = No' and M = 0) on those portions of the hinge 
lines which have transverse deflections greater than the plate thickness. 
In this case, a time dependent rectangular boundary travels outwards 
from the central line hinge toward the edge of the plate. This boundary 


always has a deflection W = H and divides the plate into two regions: 


ice 





an inner zone with W > H_ and an outer zone with W <H. The dis- 
Sipation relations (Nw - M) om must be evaluated in the inner region 
with N = No', M =0 while in the outer region the evaluation which 
was used earlier remains unchanged. Figure 4 indicates the two 


regions: 





figure 4 
If it is assumed that W/H exceeds one when t ye , then the first 


stage of motion given by equations (9a, b, and c) remains unchanged, 
while the second stage remains valid until t = tp , when W/H = 1 
where ty is the response time of the second stage alone. Thus, the 
non-dimensional time Po ot2- may be found from equation (11) with 


W =H, i.e.: 
2 ag a W ee = 
Peemepe Ge Ps 5 Po FW cprwi-H=0 (13) 


The smallest non-negative root of (13) gives Po» i.e. the response time 
for the second stage. 

Now all subsequent behavior of the plate will have Wy H at the 
centrally located plastic hinge. It is therefore necessary to consider a 
third stage of motion which is governed by the dissipation relations 
employing N = No', M =0 where necessary, as discussed earlier. It 
is quite straightforward to snow that in this case the following equation 


is obtained: 
(see next page) 
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te rr A 
We We 2 te. WW t) = 7 ue P =Q (14) 
3 1 3 Ww 4 W 
5 5 
Again the simplification p = 1 is employed so that (14) becomes: 
WW (Go Cae a 27, = 0 (14a) 
SNR 3. 3 Aus 1° 3 a is 
+ = 
feet tS T+t. pz 0 
z,=M2_, Z,:3_ , Z,- 34 pelle 
a ANG at US 
eC? = 30) : 
Pe (3-0) ’ Jo = b tang 
TW, = 2 (3 - 3 $0 +802) 
Hage) 6) 
TT . = Hb (cote +tand) 
3 Gas becota) 
- 1] 
, § (2 sing “/p 
T, = Cy QSOS ee | 
H (1 +b cot ) 
T= C,bH [cot — (tan ) a + tand |] 
: 3 : (Gist baeota> >) 
A series solution to equation (14a) with the initial conditions Ww, =H, 
Ww, = W, at te = 0, where tg is refering to third stage alone, gives: 
a a, er 
= +O + @ 15 
We H+W,t, 9 be gta (15a) 


~20- 





where: , 


il 2 2a 
a | 7 + + 4 
2 FAL ‘i A» a 43 es ) We} 
oe =. s=jzZ HW tL HW 2.a WW sa As aa fe) 
3 Sia! 1 Z 5 Z Ze 4 5 
d i, = Wi ce plod waeancy, care 
an 3 7 9 5 be 3 ts (1555) 
when te = T, ; Ws =0, Ww. = Ww, 


Substituting for T, in equation (15b) and retaining its smallest 


1 
Mon-nepative root,l,; is obtained, i.e. the response time of the third 
stage alone. The value of the final and permanent displacement is 


given then by equation (15a) with tg = T,. 


ae . ; 
> =1+ Wo ly a. SO) + Gs Ty? (16) 
H H Jal 


where Wo is obtained from equation (13) by differentiation: 


ae De SaaS 
Wie 4 2 OL + TE + 
5 ; 2 ot Po 3a, Po 4Q,T p (16a) 
An alternative approximate solution to equation (14a), for the third stage 


of motion, which is valid for any value of p, but for only a short range 


of values of the aspect ratio b, can be obtained as follows: = 


3 
Call wy, = s 
f = ' = 
We, / Zone Z, 3 HZ, 
Wb WE 
71 ye 7 ae 
2 H ’ Al H 





Z' vee 

1 2 
aoe 2 
3 4 


Then equation (14a) becomes 


ee 7 wl ag 
+Z! Xx. + lp) + + = 
a Z',W, ( 1 Wy Dy 4 (XK, W, p)=0 (17) 
W 
4 

By inspection the constants x and Xo, which both depend on the 
aspect ratio b, are approximately equal for 0 fb Go 3, their 
ratio varying between 0.98 and 1.16. 

Thus, by the substitution X; =X%, =&, equation (17) is further 


Z 
transformed to: 


+ ee NY & Sek 4 =0 Chay) 
(XK +w,"/p) 4 


Employing the change of variables 


W,=W : Wa and integrating: 


. + 
d W, 


alee ee we oleae 
wit, aw _ fw, pdWyal(W, “'p dW, 


: 1-9; 
ae ne dW, )} + 


otf wer 4h aw -b ay PO ay om 
4 





mea toe ae WC = 0 | (17b) 


For mild steel, when p = 5, equation (17b) yields: 


(see next page) 
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1 x ee 2 a, KD 68K ty 8/5 
2 no TG 1a) ors 3 Re ae 
2 32 Vy ae 
ANS ile ea ; 3) ae Comex < * eee 8/s, te 
+ (5 7 WwW, (— 5 a, 
Bea ie Tec Wee he bp 
= + o 
ae pe 2 4 Meare 
The constant C is evaluated with the conditions: 
H H 
Then the conditions 
e Wf 
att=T, Ww, = 0, Ww, “Fa give: 
Ww 
a3 ( me - + Z' log Pa ae = le ame . Te 
2 H 4 e H 2 H 
34 = 53 = 
ro -_— Z 8 
Moe 85, 2,5 . 8A 5 Sax), 2,5, 
4 9 H 8 3 H 
SW = ial =i 
a9 2 aon W 7 -— 2 W 6 /,- 
5 Nee 8 2 Os 5 f& 5 2 5 
Be f= CX —— +(— e508 OK oo + 
(“> 7 ) ( H ) (“ 5 ) ¢ H ) 
W Z 
ea ice”) ( : ee - (18) 


The final permanent deflection is given by the root of transcendental 
equation (18) and can be expected to be within the interval: 


(see next page) 
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DISCUSSION 


The approximate method presented in this pener has been used to 
predict the permanent transverse deflections (Wf/H) presented in 
Figures (5-9) for fully clamped, rate-sensitive rectangular plates 
which are subjected to uniformly distributed impulsive loading, i.e. 
Pen pic = Vo. 

It is evident from figure (5) and tables (2-6, 9) that the estimates 
which were obtained using equations (11) or (16) and rectangular plates 
with b = 0.593 agree reasonably well with the corresponding experimen- 
tal values. 

It is apparent from figures (6- 8) that, for a given value of 
impulse, the permanent transverse deflection of a plate is essentially 
independent of the magnitude of Hy when N is larger than 90. It is also 
noticed that for increasing values of n the deflections slightly decrease, 
whereas in the non-rate- sensitive analysis of Ref. [1] they are seen to 
slightly increase. The probable explanation for this difference in be- 
havior could be that in the rate- sensitive case with increasing n the 
rate- sensitive effects become more pronounced the net result being the 
(favorable) reduction in Wf/H. However, since for large values of icy 
the deflections are essentially independent of nm, as discussed above, 
for impulsive loading, when n-—»+ oo , this difference in behavior 
becomes insignificant. 

From figures (9, 10) and tables (9-11) it is observed that the 
value of the permanent deflection is also affected by the value of the 
hinge width and the value of the material property D. It should be pointed 
out here that the value of D is merely a value used in equation (2) to 
Sescribe the 6 —€ curve for the material. The proper value of D for 
mild steel is 40 se when p = 5. The mathematical difficulties that 
led to the linearization of equation (2) by setting p = 1 naturally have an 


effect on the value of D as well, thus affecting the final results. However, 


SOG. 





it is evident from figure ( ) that for values of D greater than 300 sec” 


the variation in Wf/H is rather small so that if D = 300 oe is used 
with the present analysis results good enough for engineering predic- 


tions can be obtained. Similarly, the optimum value of the hinge width, 


which is unknown at the outset, may be seen from figure ( ) to be 1 = 5H. 


(These observations are made for the particular case of mild steel). 

All the theoretical results presented in this paper were obtained 
by assuming that angular changes across the plastic hinges were suf- 
ficiently small to permit tan@ tobe replaced by 9 rad. This simpli- 
fication provides a good approximation when Wf/H is not too large. 
Moreover, only the few first terms were retained in the series solutions 
of the differential equations encountered in the analysis. It would be 
necessary to consider additional terms for large Wf/H ratios and for 
cases in which the dynamic to static pressure ratio (1) is small. The 


yield condition derived in ref [4] and presented here as equation (2) is 


1 


an approximation to the exact one, though the effect of this approximation 


is not considered to be important. The linearization of equation (2) by 
setting p = 1 is considered to be of Some importance since then a 
highly non-linear curve is replaced by a straight line. However, given 
the fact that the most important factor in dynamic loading of plates is 
the influence of finite deflections, which are ignored in bending- only 
methods, this linearization seems to be acceptable,at least tentatively, 
Since the results are reasonably close to experimental values and at 
the same time a good insight into the rate-sensitive behavior is gained. 
The approximate procedure presented in this paper could be 
used to examine the influence of large transverse dynamic loads on the 
behavior of beams and plates which have any shape and support condi- 
tions, provided the appropriate plastic hinge pattern is postulated. 
Finally it should be remarked that the predictions of this analysis 


are thought to be valid when the external dynamic energy applied toa 


-26- 





beam or plate is appreciably larger than the amount of energy which 


could be absorbed in an elastic manner. 


CONCLUSIONS - RECOMMENDATIONS 


An approximate theoretical investigation is developed herein in 
order to estimate the permanent deflections of rectangular plates and 
beams subjected to large dynamic loads. The influence of geometry 

changes and of strain-rate sensitivity is retained in the analysis but 
Hoeclastic effects are considered. The particular case of a fully 
clamped rectangular plate acted on by a uniformly distributed pressure 
pulse is studied in detail. Tables and graphs are provided for the case 
of fully clamped plates and beams made from mild steel and loaded 
impulsively. It is observed that good agreement between the theoretical 
predictions and experimental results has been obtained for a rectangular 
plate of aspect ratio b = 0.593. Improved theoretical predictions are 
expected to be obtained if equation (11) is solved for the correct value 
of p. Though equation (18) is an exact solution to equation (17), it is 
valid only for a short range of aspect ratios (0 - 0.3) and even in that 
range it contains some approximations. Moreover, equation (17) 
corresponds to the third stage of motion only so that the exact solution 
(18) should not be used unless equation (11) has been solved correctly 
as well. Exact or numerical solutions to equations (11) and (16) would 


be a useful complement to this paper. 
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(11)or(16) Ref [2] 
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Table (2). b=0.593, 1=5H, D=300sec , n = 200. 
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Equation Experimental, 
(11)or(16) Ref [2]. 
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Table (3). Theoretical vs. experimental results. b = 0.593, D = 460 ee 
ifaezeod, 1 = 200. 
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Equation Experimental 
(11)or(16) Ref [2] 


H-in Vo- ft/sec Wf/H Wf/H 
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Equation Experimental 
(1l)or(16) Ref [2] 
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Equation Experimental 
(11l)or(16) Ref [2] 


D-sec : H-in Vo-ft/sec Wwt/H Wf/H 
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Table (9). Influence of D on predictions based on equations (11) or 


(i6ys 1 = 2.5H, y = 200. 
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Table (10. Influence of 1/H on predictions based on equations (11) or 


(16). b= .0593, m = 200, D = 470 ae ee . 0643, = 417.16. 
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Table (1). Influence of hinge width on the predictions based on 


equations (11) and (16). b =0.593, D=500sec . m= 200. 
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